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Abstract 

A Fermi resonance-algebraic model is proposed for molecular vibrations, where a 
U(2) algebra is used for describing the vibrations of each bond, and Fermi resonances 
between stretching and bending modes are taken into account. The model for a bent 
molecule XY2 and a molecule XY3 is successfully applied to fit the recently observed 
vibrational spectrum of the water molecule and arsine (AsHa), respectively, and 
results are compared with those of other models. Calculations show that algebraic 
approaches can be used as an effective method for describing molecular vibrations 
with small standard deviations. 

Keywords U(2) algebra Vibrational spectra Fermi resonance 



1 



Algebraic methods have been appUed to molecular systems for a number of years, and 
a recent review is available. U(4) and U(2) algebraic model have mostly been used so 
far for a description of rovibrational spectra of molecules. U(4) model has two advantages 
that rotations and vibrations are treated simultaneously, and that Fermi resonances are 
described by the nondiagonal matrix elements of Majorana operators. But it becomes 
quite complicated when the number of atoms in a molecule is larger than 4. U(2) model is 
particularly well suited for dealing with the stretching vibrations of polyatomic molecules, 
and currently under further development 

Recently, an important development on these methods is to incorporate bending vibra- 
tional modes of polyatomic molecules. Sample molecules such as acetylene as well as 
X3 and XY4 [^''''1 molecules have been studied via algebraic methods. These approaches 
suffer from the neglect of the Fermi resonance interactions between stretching and bending 
modes. Due to their importance as a mechanism for intramolecular energy transfer and 
for descriptions of excited stretching and bending vibrations, Fermi resonances have been 
taken into account in simple Fermi resonance-local mode models for bent triatomic I^' and 
pyramidal XH3 systems and in boson-reahzation models for bent triatomic 1^°! and XY4 
[^^1 molecules. In this paper we will introduce a Fermi resonance-algebraic model, where a 
U(2) algebra is used for describing the vibration of each bond. Our method is addressed 
for a bent molecule XY2 and a molecule XY3, where Fermi resonances between the stretch 
and the bend are considered. In a limit, the model corresponds to the boson-realization 
model. Vibrational spectra of two sample molecules H2O and AsHa are used to test our 
model. This method can be extended for describing the vibrations of molecules with other 
symmetry. The corresponding results will be discussed elsewhere. 



1. Hamiltonian 

Let n be the total degrees of freedom of vibrations and m the degrees of freedom of 
stretching vibrations in a polyatomic molecule. We introduce n U(2) algebras to describe 
vibrations: Ui(2) (1 < i < m) for stretching modes and U^(2) {1+m < n <n) for bending 
modes. Each Uq(2) (1 < q; < n) is generated by the operators {TVq,, J+^aJ-,a: Jo,a}, 
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satisfying the following commutation relations: 

[Jo,a, J±,l3\ — ±(^a/3-^±,a, [J+,a, J-,l3\ — '^3al3Jo,a, 

[A^a, io,/3] = 0, [A^a, 4/3] no- 

where Na is related with the Casimir operator of U(2): 

2Jo% + J+,aJ-,a + J-,aJ+,a = Na{Nj2 + 1). 

The local basis states for each bond are written as IA^^q,, Va), where is the eigenvalue 
of Na and Va the number of quanta on the ath bond. The action of J±^ce on the local 
states is given by 



J+^a\Na,Va) = VaiN^ - Va + I) {N^, - I) , 
J_,Q \Na,Va) = \/{Va + l){Na- Va) \Na,Va + l). 

Those Na of equivalent bonds are equal to each other: Ni = Ng, and A^^ — N^, where 
footnotes s and b refer to the stretching and the bending vibrations, respectively. 

Due to Fermi resonances in the considered system, we will take V — Vs + Vf, /2 as a, 
conservative quantity, where Vg and Vj, is the total number of phonons of stretching and 
bending modes, respectively. Define the scale transformations for convenience: 

Oa = J+,a /\fN^, ai = J-,a / ^j^a- (1-3) 

In the limit of A''^ ^ oo and A^^ — >• cxd, aa(aj^) reduce to bosonic operators with usual 
boson commutation relations '^l 

In the following, we will present the vibrational model in terms of Eq. (1.3) for two 
molecular systems with different symmetries. 

1.1 The model for vibrations of a bent molecule XY2 

For a bent triatomic molecule, we have n — 2) and m — 2. The Hamiltonian, where 
the interactions are taken up to the fifth orders, is given as follows: 

2 

if = ^ a\ai{us + Xsttia]) + 403(0;;, + xi,a^a\) + \i{a\a2 + a\ai) 

i=l 

+ X2{a\aia2a2) + X3{a\a\a2a2 + H.c.) + X4{a\a2 + a2ai)(a|ai + ^02) 

+ A5(a|a2 + 0201)0303 + A6(o|oi + 0202)0303 + AyKoj + 02)0303 + H.c.} ^^'^"^ 

+ A8{o|o2(oi + 02)0303 + if.C.} + AgKojojoi + 020202)0303 + i/.C.} 

+ Aio{(aiaia2 + 020301)0303 + H.c.} + Aii{(ai + 03)03030303 + H.c.}, 
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where uj and x arc used for comparison with the constants in the Morse potential Aj 
(1 < i < 11) are the couphng constants. The term with A7 describes the Fermi resonance 
between the stretching and the bending vibrations. The last four terms in Eq. (1.4) are 
the fifth order interactions caused by Fermi resonance. In the hmit of Ng — > 00 and Nb 
— > 00, Eq. (1.4) corresponds to the boson- reahzation model I-^^l. 

1.2 The model for vibrations of a molecule XY3 

For a C-iv symmetric molecule, we have n = 6 and m = 3. From the standard 
method of group theory, it is not difficult to obtain the character table, the representation 
matrices of the generator, and the Clebsch-Gordan coefficients of the Cs^, group. From 
that knowledge, Hamiltonian where the interactions are taken up to the fourth orders are 
expressed in terms of operators Eq. (1.3): 

3 6 
i=l M=4 iy^j 

+ V2 Y aj^aj, + r]3 {a\a4a4^ + ala^a^ + alaQttQ + H.c.) 

+ 774 {01(0505 + oeOe) + 02(^404 + oeOe) + 03(0404 + 05^5) + H.c.} 
+ rjs {o|o4(o5 + oe) + 0^05(04 + Og) + 0305(04 + 05) + H.c.} 
+ r]Q (aJasOe + 020405 + 030405 + H.c.) + ^77^ o|oio]oj 

+ Vs Yl ("'hi + + + V9 Yl alai{a]ak + ojo^) 

+ ^710 Yi i^hl^j^j + ^■^■) + ^11 m («i«Jai«fe + H.c.) 

(1.5) 

where cu and x are again used for comparison with the constants in the Morse potential. r)i 
and 7)2 is the harmonic constant for the stretch and the bend, respectively. These rjp (3 < 
p < 6) are the Fermi resonance coupling constants, i]q {7 < q < 11) anharmonic constants 
of the stretching vibrations, and i]r {12 < r < 15) anharmonic constants of interactions 
between the stretch and the bend. In Eq. (1.5), we have omitted the interaction terms 
for the bending vibrations, which are similar to the terms with 77^ (7 < < 11) for the 
stretching modes, and also omitted the fourth-order interactions between the stretching 
and the bending vibrations. In the limit of — > 00 and N;, — > 00, equation (1.5) reduces 
to the extended local mode model for three equivalent stretching bonds '^^1, where the 
bending vibrations were neglected. 
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2. Applications 



2.1 Application to H2O 

Studies of the vibrations of the water molecule already exist in the literature. Its 
stretching vibrational energy levels were calculated via nonlinear quantum theory '^^1. 
Halonen and Carrington proposed a simple vibrational curvilinear internal coordinate 
hamiltonian t^l for its vibrational levels up to 18,500 cm~^. The same levels were studied 
by U(4) algebraic model '^l Its recently observed vibrational spectrum up to 22,000 cm~^ 
were analyzed by the potential-energy surface from ab initio calculation '^^1 as well as in 
the boson-realization model and the corresponding gf-deformed model Owing to much 
research on this molecule, it provides a good testing ground for different models. 

Now, we calculate the Hamiltonian matrix elements in the symmetrized bases, then 
determine the 15 parameters in Eq. (1.4) by a least-squares optimization in fitting the 
experimental data Eobs in cm~^ (their normal labels (z^i^'2K3)), taken from the compilation 
of Ref. [14]. Two boson numbers and are taken to be 48 and 65, respectively. The 
parameters obtained are given in cm~^ as follows: 

a;, = 3705.121, = -4.602, = 1599.485, a;,, = 3.913, Ai = -51.335, 
A2 = -12.746, A3 = -1.004, A4 = 3.199, A5 = -2.853, Ae = -21.956, 
A7 = 12.686, Ag = -0.422, Ag = -5.901, Aio = 2.187, An = 4.220. 

Table 1 

The observed and calculated vibrational levels arc presented in Table 1 where our 
results are also compared with those calculated in the boson- realization model ^^'^1. The 
AEcai are the differences between the observed and the calculated values. The standard 
deviation (SD) of this fit is 6.724 cm~^, while that of the boson- reahzation model and 
the corresponding g-deformed model is 8.198 cm~^ and 7.157 cm~^ 1^°!, respectively. The 
present model applied to H2O gives a slight improvement in the fit to its vibrational 
energy levels. 

2.2 Application to ASH3 

To our knowledge, vibrations of a molecule XY3 have not been studied via algebraic 
methods yet. As an example, we now apply Eq. (1.5) to the vibrational spectrum of arsine 
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(AsHa). The calculation of energy levels will be greatly simplified if the symmetrized bases 
are used. For the considered system, these bases can be obtained by the method of Ref. 
[7] . This method is particularly useful when describing vibrations of polyatomic molecules 
and for high overtones. 

Recently observed data (-Bobs) AsB.^ with their normal labels are taken 

from the compilation of Refs. [9,15]. Since there are 18 experimental data for this molecule 
it is not possible to determine all coupling constants in Eq. (1.5) by a least-squares 
optimization. It is natural to choose the following parameters in Eq. (1.5) nonvanishing: 
cUs, Xg, cUb, Xb, r)i, r)2, 773, 774, 775, and 776- In other words, we neglect the terms with 77^ 
(7 < i < 11). We have found out that other choice does not provide improvement in our 
fits. For comparison, we have made the first fit (Fit a) where the model is in the limit of 
Ng ^ 00 and N^, ^ 00. The obtained parameters and the calculated differences (AE^^i) 
with the corresponding observed data (-Eot^) in this fit are listed in Table 2 and 3. In 
another fit (Fit b) we take the boson numbers A^^ and Nb to be 62 and 42, respectively. 
The results in Fit b are also fisted in the related Tables. 

Table 2 
Table 3 

The SD in Fit a and b is 10.201 cm~^ and 5.827 cm~^, respectively. The calculated 
results for this molecule in Fit b can be compared with those by Lukka et al '^1 . However, 
their vibrational Hamiltonian was expressed in terms of curvilinear internal coordinates, 
and the Morse oscillator basis functions and the results of ab initio calculation were used 
there. 

3. Conclusion and discussion 

We have proposed Fermi resonance- algebraic models for the vibrations of a bent 
molecule XY2 and a molecule XY3 by U(2) algebras for the description of both stretching 
and bending modes, where Fermi resonances between the stretch and the bend are consid- 
ered. In the limit, the models reduce to the recently presented boson-realization models. 
As an example, they have been applied to vibrational spectra of the water molecule and 
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arsine (AsHa). To some extend, they have provided improvement in the standard devia- 
tions in the corresponding boson-reahzation models for the same observed energy levels. 
It needs more sample calculations to judge which approach is better. 

The present calculations and others [^^'^''^^'^^l demonstrate that an algebraic approach 
can be used as an effective method for describing molecular vibrations with good precision. 
In this approach the eigenvalues and the related wave functions are obtained through 
matrix diagonalization. Hence, the required computing time is short, and some other 
physical properties, such as transition intensities t^'^^J, are allowed to calculate. It provides 
the possibility of using an algebraic model as a numerically efficient, phenomenological 
theory for studying molecular spectra, especially when no ab initio calculation is available. 
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TABLE 1. Observed and calculated vibrational energy levels of H2O ( cm ^) 





fl4l 
TP'- ' 






(1/11/2^3) 








(010) 


1594.7498 


-2.668 


-4.733 


(001) 


3755.93 


0.170 


-0.522 


(020) 


3151.63 


0.016 


-4.136 


(Oil) 


5331.269 


-5.434 


-5.562 


(100) 


3657.053 


3.549 


1.997 


(021) 


6871.51 


-3.398 


-4.325 


(030) 


4666.793 


5.464 


0.110 


(101) 


7249.81 


2.611 


0.409 


(110) 


5234.977 


1.684 


0.694 


(031) 


8373.853 


4.051 


2.443 


(040) 


6134.03 


9.088 


4.828 


(111) 


8807 


-1.858 


-0.673 


(120) 


6775.1 


1.817 


-0.591 


(041) 


9833.584 


13.643 


12.564 


(200) 


7201.54 


1.302 


0.483 


(121) 


10328.731 


-0.616 


-0.006 


(002) 


7445.07 


3.063 


1.453 


201) 


10613.355 


0.766 


-0.809 


(050) 


7552 


11.254 


10.415 


(003) 


11032.406 


5.575 


3.577 


(130) 


8273.976 


0.594 


-1.912 


(131) 


11813.19 


2.636 


1.919 


(210) 


8761.582 


0.248 


1.828 


(211) 


12151.26 


-3.666 


0.371 


(060) 


8890 


-17.089 


-12.506 


(013) 


12565 


-9.461 


-6.151 


(012) 


9000.136 


-5.600 


-4.650 


(141) 


13256 


4.292 


4.464 


(220) 


10284.367 


-0.405 


-0.714 


(221) 


13652.656 


-3.450 


-0.116 


(022) 


10524.3 


-0.950 


-0.991 


(301) 


13830.938 


0.384 


-1.448 


(300) 


10599.686 


-3.117 


-3.002 


(023) 


14066.194 


-9.806 


-7.824 


(102 


10868.876 


5.583 


1.810 


(103) 


14318.813 


6.331 


1.985 


(310) 


12139.2 


-6.313 


-1.508 


(151) 


14640 


-10.977 


-5.968 


(112) 


12407.64 


3.124 


3.469 


(231) 


15119.029 


-3.536 


-2.852 


(042) 


13448 


7.777 


4.483 


(311) 


15347.956 


-4.276 


1.217 


(320) 


13642.202 


-7.521 


-4.024 


(113) 


15832.765 


-7.239 


-1.381 


202 


13828.277 


2.052 


0.344 


321) 


16821.635 


-2.039 


2.715 


(122) 


13910.896 


6.188 


-0.130 


(203) 


16898.842 


3.234 


-1.570 


(400) 


14221.161 


-1.522 


-1.806 


(123) 


17312.539 


1.834 


-1.041 


(004) 


14536.87 


9.844 


6.393 


(401) 


17495.528 


0.756 


-0.563 


(330) 


15107 


-14.579 


-12.251 


(331) 


18265.82 


-3.654 


-2.516 


(212) 


15344.503 


12.518 


10.038 


(213) 


18393.314 


-1.400 


0.480 


(410) 


15742.795 


0.179 


4.839 


(411) 


18989.961 


-12.986 


-0.621 


(222) 


16825.23 


-16.629 


-18.041 


(421) 


19720 


-0.801 


-0.492 


(302) 


16898.4 


1.469 


-3.711 


(303) 


19781.105 


14.596 


8.509 


(420) 


17227.7 


11.514 


0.715 


(501) 


20543.137 


-7.555 


-2.973 


(500) 


17458.354 


0.824 


0.139 


(313) 


21221.828 


-6.455 


-5.430 


(104) 


17748.073 


-0.279 


-0.504 


(232) 


18320 


25.584 


29.589 


(312) 


18392.974 


-1.070 


0.979 


(412) 


21221.569 


-5.258 


-5.451 












SD 


8.198 


6.724 
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TABLE 2. Parameters (cm ^) for AsHa obtained by the least square fitting 



Fit a 2039.921 
Fit b 2066.972 



Xs Ub Xb r]i 

-80.580 959.021 64.320 -24.775 
-28.376 962.647 68.865 -18.639 



m V3 V4: V5 V6 

Fit a -21.932 -53.616 -5.996 -33.712 1.823 
Fit b -27.283 -55.482 -5.850 -29.862 11.970 



TABLE 3. Observed and calculated vibrational energy levels of AsH^ (cm ^) 













^obs 






(OIOOqO) 


906.752 


-8.425 


-1.328 


(000° l±i) 


999.225 


8.253 


9.295 


(02000°) 


1806.149 


3.950 


1.853 


(010°l±i) 


1904.115 


-9.368 


-0.427 


(00002°) 


1990.998 


-2.224 


1.602 


(000°2±2) 


2003.483 


-10.736 


-9.707 


(IOQOoO) 


2115.164 


-0.720 


0.932 


(OOI^OO) 


2126.423 


2.698 


-1.537 


(IIQOoO) 


3013 


1.704 


3.953 


(ooiH±i) 


3102 


-0.615 


3.396 


(200^00) 


4166.772 


-0.182 


3.033 


(lOl^oO) 


4167.935 


-8.389 


-2.245 


(30000°) 


6136.316 


2.401 


-4.667 


(20110°) 


6136.310 


9.419 


1.951 


(00300°) 


6276 


-3.991 


-2.321 


(102^0°) 


6295 


2.410 


2.238 


(40000°) 


8028.977 


-2.275 


-1.389 


(301^0°) 


8028.969 


1.398 


0.713 



SD 10.201 5.827 
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